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Abstract. The accurate characterization of the nuclear symmetry energy and its density 
dependence is one of the outstanding open problems in nuclear physics. A promising nuclear 
observable in order to constrain the density dependence of the symmetry energy at saturation 
is the neutron skin thickness of medium and heavy nuclei. Recently, a low-energy peak in the 
isovector dipole response of neutron-rich nuclei has been discovered that may be correlated with 
the neutron skin thickness. The existence of this correlation is currently under debate due to 
our limited experimental knowledge on the microscopic structure of such a peak. We present 
a detailed analysis of Skyrme Hartree-Fock (HE) plus random phase approximation (RPA) 
predictions for the dipole response in several neutron-rich nuclei and try to elucidate whether 
models of common use in nuclear physics confirm or dismiss its possible connection with the 
neutron skin thickness. Finally, we briefly present theoretical results for parity violating electron 
scattering on ■^°*Pb at the conditions of the PREx experiment and discuss the implications for 
the neutron skin thickness of ^"^Pb and the slope of the symmetry energy. 



1. Introduction 

Experimental studies [H [2l [3] on the low-energy isovector dipole response or Pygmy Dipole 
Resonance (PDR) in neutron rich nuclei are of crucial relevance because they determine reaction 
rates in the r-process [3| and, in addition, the PDR has been related to the neutron skin thickness 
of the studied nuclei in Refs. [5l|6l[7]. However, the existence of this correlation is currently 
under debate [8] due to our limited experimental knowledge of the microscopic structure of 
such a peak. If the underlying dynamics giving rise to the PDR is, eventually, confirmed to 
be strongly correlated to the formation of a neutron skin in neutron-rich nuclei, the density 
dependence of the nuclear symmetry energy at saturation will be immediately constrained [Sl\- 
This is of broad interest since the symmetry energy and its density dependence impact on a 
variety of physical systems, such as the composition and structure of the crust in a neutron star 
[lOl [11] , the neutron skin thickness of a heavy nucleus O [12l [13] , atomic-parity violation [H] 
and heavy ion collisions [15j. 

For all these reasons, we have performed in Ref. |16j a detailed analysis of Skyrme Hartree- 
Fock (HF) plus random phase approximation (RPA) predictions for the isovector and isoscalar 



dipole response in ^^Ni, ^^^Sn and ^'^'^Pb nuclei — representative of different mass regions — in 
order to elucidate the nature and possible connection of the PDR with the slope of the symmetry 
energy. The strategy adopted to understand if such a connection may exist, is very simple. 
Giant resonances are collective excitations of atomic nuclei that have allowed, in the past, to 
determine some nuclear saturation properties such as the nuclear incompressibility (strongly 
related to the Giant Monopole Resonance [E]), the effective nucleon mass (which affects the 
Giant Quadrupole Resonance [IE]), or the nuclear symmetry energy at sub-saturation density 
(which act as a restoring force in the Giant Dipole Resonance [12]). Hence, collective phenomena 
inform us about general properties of the nuclear effective interaction and all realistic models 
should predict such a collectivity even though some differences on details may appear. The 
additional fact that the energy weighted sum rule for the PDR has been correlated with the L 
parameter [6], defined as L = 3po where Csym(p) is the symmetry energy, p the baryon 

density and po the nuclear saturation density, has motivated our selection of studied interactions. 
Specifically, we use three Skyrme interactions widely used for nuclear structure calculations and 
that differ in their predictions of the L parameter to study the collectivity displayed by the 
RPA states giving rise to the pygmy dipole strength, or RPA-pygmy state, as well as its possible 
relation with the neutron skin thickness. Based on our experience, we present here the common 
features found in different mass regions by using 208p]-, ^ representative example. 

Currently, the PREx collaboration [20j aims to determine the neutron radius of 208p]-, ^j^hin 
a 1% error by parity violating electron scattering (PVES) [21]. Such a measurement is very 
important for three basic reasons. First, it measures the neutron distribution in a heavy nucleus 
free from most of the strong interaction uncertainties. Second, it paves the way for further 
measurements of neutron densities by PVES |23l I24j . Third, it may allow one to derive a 
significant constraint on the L parameter [22l[25] and, therefore, it may help in constraining the 
isovector channel of the nuclear effective interaction. Here, we will shortly discuss mean-field 
model predictions for PVES at the kinematics of PREx j22] . 

The work is organized as follows. In Sec. [2] we briefly present the basic formalism employed 
in our analysis. For further information we address the reader to Refs. I26j . In Sec. [3l we 
show the main results of our works \\.&\ 122]. Finally, our conclusions are laid in Sec. |H 

2. Formalism 

2.1. Random Phase Approximation 

The RPA method is well-known from textbooks |26] . In short, once the HF equations are 
solved self-consistently for the given Hamiltonian, we build accordingly the residual interaction 
— considered to be a small perturbation of the HF mean-field potential — and, then, we solve the 
RPA coupled equations by means of the matrix formulation. The important quantities to define 
for our study are the following. The main one is the reduced transition strength or probability. 
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where Aph{EJ,0 — > is the reduced amplitude, |0) is the RPA ground state, {u) is a generic 
RPA excited state, J is the angular momentum carried by the operator, Fjm, meant to modelize 
the experimental probe and is the reduced matrix element of such an operator between 

a hole (/i) state (occupied state) and a particle (p) state (unoccupied state). The sum of all 
ph states that contribute to an RPA transition is weighted by the X'^ and Y'^ amplitudes, 
eigenvectors of the RPA secular matrix L26j. For further details we also refer to 



2.2. Parity Violating Elastic Electron Scattering 

Parity violating electron-nucleus scattering (PVES) probes neutrons in a nucleus via the 
electroweak interaction [281 ES] ■ Electrons interact with the protons and neutrons of the nucleus 
by exchanging a photon or a boson. The former mainly couples to protons while the latter 
basically couples to neutrons. In the experiment, one measures the parity- violating asymmetry, 
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where da±/dn is the elastic electron- nucleus differential cross section for incident electrons with 
positive and negative helicity states. For a realistic calculation of the parity violating asymmetry, 
we solve the Dirac equation via the distorted wave Born approximation (DWBA) [3D] where the 
main input are the electric and weak charge distributions of the studied target 



3. Results 

3.1. Low-energy dipole response of'^^^Pb 

In our recent work we have studied in detail the dipole response of ^®Ni, ^^^Sn and 208p^ 
by means of the formalism explained in Section [2.11 and by using different Skyrme interactions. 
Here we will present the case of 208pb as a representative example of the common trends found 
in different mass regions. 
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Figure 1. Strength function calculated by convoluting the corresponding reduced transition 
probability [Eq. ([T|)] with a Lorenzian of 1 MeV width for the isovector (left) and isoscalar 
(right) dipole response of 208p|-, ^ function of the excitation energy [IB]. In both figures 
the predictions of SGII, SkI3 and SLy5 are depicted. Black arrows indicate the experimental 
centroid energies for the PDR (E = 7.37 MeV) [I] and for the IVGDR (E = 13.43 MeV) [S^. 

In Fig. [1] we show the isovector (IV: left panel) and isoscalar (IS: right panel) averaged strength 
functions for the dipole in ^"^^Pb predicted by the Skyrme models SGII [31j with an L = 37.6 
MeV, SLy5 [32] with an L = 48.3 MeV and SkI3 [33] with an L = 100.5 MeV as a function of the 
excitation energy. Experimental data for the centroid energies of the PDR and the IV Giant 
Dipole Resonance (GDR) |34j are also depicted (black arrows). The theoretical predictions lie 
within the experimental error and, therefore, the Skyrme-HF plus RPA approach may constitute 
a good starting point for a detailed analysis of the microscopic structure of the PDR in those 
nuclei. An interesting feature that can be also observed in different mass regions [16] is the ratio 
between the low- and high-energy strength exhausted by the IS peaks: of the same order; and 



by the IV peaks: one order of magnitude. This aheady indicates the importance of the probe 
used to excite a certain RPA state and reveals also its nature. In other words, the IS dipole 
operator excite more efficiently the RPA-pygmy state than the IV dipole operator. This means 
that a perfect probe for exciting such a state would be mostly isoscalar. Moreover, we see from 
both figures that as the predicted value for the L parameter increases — in going from SGII to 
SLy5 and, finally, to SkI3, the low-energy peak is shifted to larger excitations energies and to 
larger strengths, in perfect agreement with Ref. 
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Figure 2. Neutron (black) and proton (red) ph contributions to the isovector (left) and 
isoscalar (right) reduced amplitude [Eq. ([T])] corresponding to the RPA-pygmy state of 208p|-, 
(RPA excitation energy is indicated) predicted by the different models as a function of the ph 
excitation energy pjB]. The single particle levels involved in the most important ph transitions 
are indicated. 

In Fig. [21 we show the ph contributions to the IV (left) and IS (right) reduced amplitudes 
[see Eq. ([1])] corresponding to the different RPA-pygmy states of 208p]-, predicted by the different 
models as a function of the ph excitation energy. Two results — common also to other mass 
regions — arise from the analysis of this quantity. First, while the isoscalar reduced amplitude 
{Apfi{El, IS)) is formed by several neutron ph transitions adding coherently (same sign) and 
only few states are contributing destructively, the isovector reduced amplitude {Aph,{El, IV)) 
is formed by few neutron and proton ph transitions contributing with different signs, i.e., 
non-coherently. Second, the most relevant ph transitions in the isoscalar dipole response of 
208pj^ correspond to excitations of the outermost neutrons that, as a consequence, dominate the 
dynamics in the low-energy region (see [E]). For the case of the isovector dipole response the 
situation is model dependent. However, the most relevant neutron ph contributions are larger 
in number than the proton ones and basically due to the outermost neutrons. Therefore, albeit 
some differences will be present, one may expect a relation between the neutron excess and the 
IS and IV dipole responses in neutron rich nuclei. 

From Fig. [H we have seen that the low-energy IS and IV dipole responses of ^"^Pb display, to 
different degree, a sizeable low-energy peak in the corresponding strength functions (see section 
III.B of Ref. |I6]). This is actually one of the main characteristics one asks to a collective state. 
The second one is that collective phenomena (or resonances) should display coherence between 
the contributions of several ph transitions to the reduced amplitude. Therefore, by looking at 
Fig. [21 we can state that while all models support a clear collective character of the low-energy 
peak in the IS dipole response in ^"^Pb (also true in other mass regions), the collectivity of its 
IV counterpart is model dependent. 



3.2. Parity violating electron scattering on ^°Ph 

The Lead Radius Experiment (PREx) at the Jefferson Laboratory has recently reported first 
results for PVES on ^ospb [20j. In this first run statistics were not sufficient in order to achieve 
the desired accuracy ^20j. A second run of PREx has been approved and it is intended to be 
performed in the future pOj [23], 

As stated in the introduction, the relevance of this measurement has motivated our study of 
the parity violating asymmetry at the PREx kinematics |21|, [22] . We have applied the formalism 
mentioned in Section \2.'1\ which is presented in more detail in Refs. [211 [221 [211 [32] • In Fig- [31 
we display the linear correlation between the parity violating asymmetry and the neutron skin 
thickness of 208p]-, (Ar„p = (r^)^^^ — {t'^)^^'^) as predicted by more than forty mean-field models 
of very different nature: from non-relativistic such as Skyrme or Gogny models to relativistic 
such as non-linear Walecka or density dependent meson-exchange and point-couling models. All 
of them accurately reproduce the charge radius of 208p]-, ^£qj. f^pf^j-^er details see Ref. [22] )• This 
high linear correlation allows one to accurately extract the value of the neutron skin thickness 
of 208pb without assuming a particular shape for the nucleon spatial distributions. In the right 
panel of Fig. [3] we show the well-known correlation between the neutron skin thickness in 208p|j 
and the slope of the symmetry energy predicted by the considered mean-field models. Therefore, 
PVES can supply new constraints on the value of the L parameter. 




Figure 3. Left panel: parity violating asymmetry (DWBA) in 208p]-, £qj. ^ Qg Qgy electrons at 
5° scattering angle as a function of its neutron skin thickness [22] • MF results (black circles) 
— references to all these interactions can be found in [22] — and calculations form the neutron 
densities deduced from experiment [a] |35j, [b][36j and [c] [S^ (red squares). Right panel: 
Neutron skin thickness of ^"^Pb as a function of L as predicted by the same MF interactions 
shown in the left panel [22] . 



4. Conclusions 

We have exemplified some of the common features found in the dipole response of different nuclei 
[16] by showing our results for the case of ^'^^Pb. In particular, the low-energy peak in both the 
isoscalar and isovector dipole responses is shifted to larger excitation energies and display larger 
values of the strength function as the value of the L parameter increases. We have also seen that 
the RPA-pygmy state can be more efficiently excited by an isoscalar probe than by an isovector 
probe. This indicates the dominant isoscalar character of such a state. We have demonstrated 
that the low-energy peak in the isoscalar dipole response is a collective mode — it is formed by 
several ph contributions adding coherently and giving a contribution to the reduced transition 
strength comparable to that of the isoscalar giant dipole resonance. This is found to be opposite 



to what happens to its isovector counterpart where its cohectivity depends on the interaction. 
Finahy, the dynamics in the low-energy region in the isoscalar dipole response of neutron-rich 
nuclei is clearly dominated by the outermost neutrons — those that form a neutron skin. 

Our analysis of PVES applied to the conditions of the PREx experiment predicts a high- 
quality correlation between the parity violating asymmetry and the neutron skin thickness of 
^''^Pb [22]. The results suggest that one will be able to extract significant constraints on the 
slope of the nuclear symmetry energy at saturation if the statistics of PREx are improved. 

In conclusion, the recent experimental and theoretical studies of the ground-state and 
excitation properties of neutron-rich nuclei aim to complement each other and are paving the 
way for a better knowledge of the isovector channel of the nuclear effective interaction. 
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